Renormalization of f-levels away from the Fermi energy in electron excitation 
spectroscopies: Density functional results of Nd2_2;Cea;Cu04 
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Relaxation energies for photoemission where an occupied electronic state is excited and for inverse 
photoemission where an empty state is filled are calculated within the density functional theory 
with application to Nd2-a;CexCu04. The associated relaxation energies are obtained by computing 
differences in total energies between the ground state and an excited state in which one hole or one 
electron is added into the system. The relaxation energies of f-electrons are found to be of the order 
of several eV's, indicating that f-bands will appear substantially away from the Fermi energy (Ef) 
in their spectroscopic images, even if these bands lie near Ep. Similar shifts are obtained for the 
Gd-f states in Gd2Cu04. Our analysis explains why it would be difficult to observe f electrons at 
the Ef even in the absence of strong electronic correlations. 

PACS numbers: 71.15.Qe,79.60.-i,74.72.Ek 



I. INTRODUCTION 

Compounds containing rare-earth (RE) or actinide el- 
ements display many intriguing solid state phenomena 
such as heavy fermion behavior and high-temperature 
superconductivity. When partially filled f-orbitals are in- 
volved, the ground state predicted by the density func- 
tional theory (DFT) clearly places the f-electrons in nar- 
row bands piled at the Fermi energy Ep , interacting only 
weakly with other electrons. In sharp contrast however 
signatures of f-bands are often found in spectroscopic 
measurements not at Ep, as band theory predicts, but 
several eVs above or below the Ep depending on the na- 
ture of the spectroscopy [IHll. Here we show how this 
dilemma can be resolved and how this seemingly contra- 
dictory behavior of f-electrons in solid-state systems can 
be modeled within the framework of the DFT. In essence, 
we have carried out DFT calculations constrained to sim- 
ulate the process of electron excitation. In this way, we 
adduce that when f-bands are excited, their spectroscopic 
image will generally avoid Ep even when these bands lie 
at the Ep. 

We consider NCCO as an exemplar complex system, 
which is interesting not only because it is a high- Tc su- 
perconductor Q , but also because it contains Nd-f bands 
co-existing with the broader Cu-0 bands [Bj . This allows 
us to delineate relative differences in the way localized 
bands get excited in comparison to the itinerant bands. 
Our computations are based on the final state rule and 
the A-self-consistent field (ASCF) method, which have 
been invoked previously to investigate x-ray absorption 
and emission I6|-l8| , and x-ray photoemission spectroscopy 
(XPS) of core levels @. The final state rule postulates 
that the electronic structure of the excited state can be 
obtained by using the potential of the final state in which 
the hole is present in the valence band, i.e. the hole is 
long-lived and the electronic system relaxes before the 
hole recombines with an electron. In standard computa- 
tions, the excitation energies are approximated by band 



energies, but in the ASCF method the energies are sig- 
nificantly more realistic because they are obtained by 
computing the total energy difference between the un- 
perturbed ground state and a relaxed cell calculation for 
the final state. Our goal is to capture the fundamen- 
tal mechanism responsible for why f-electron excitations 
avoid Ep, without consideration of probe dependent ma- 
trix elements [Tol - [l5| . For this purpose, we introduce 
a method beyond the many body perturbation theory 
for calculating the energy cost of localized f-state 
relaxations. 



II. METHOD OF CALCULATION 

Filled states can be probed via x-ray emission or pho- 
toemission spectroscopy, and empty states through x-ray 
absorption or inverse photoemission or bremsstrahlung 
isochromat spectroscopy, issues of surface sensitivity of 
photoemission notwithstanding. Formally, the energy for 
exciting an electron requires us to compute the difference 
between the total energy Eo{N) of the TV-particle ground 
state and the energy En{N — 1) of an excited state of 
{N — 1) particles containing a hole in the n*'' level [20| . 
The relaxed excitation energy e„ of the n*^ filled level 
then is 



E,,iN -l)-EoiN). 



(1) 



The energy contribution from the photon, huj, can be 
subtracted off, since it is not an interesting part of the re- 
laxation energy. Similarly, for the inverse process, where 
we add an electron into an empty state of the A^-particle 
system, we need to evaluate the total energy En{N + 1) 
of an excited state of the (A^-|- 1) particle system, so that 
the excitation energy e„ of the n*'' unfilled state is 



E„{N + l)-EoiN). 



(2) 



In the absence of relaxation, e„ = where e° is the 



energy of the n*'* Kohn-Sham orbital [2l| . The correction 



2 



to this Kohn-Sham energy is the relaxation energy 



(3) 



associated with the excitation process. The connection 
of groundstate Kohn-Sham equations to one particle en- 
ergies has been studied by Bauer [11]. Other authors 
[Ta . [23| have shown that the exchange-correlation po- 
tential Vxc for calculating ground state properties is the 
best local approximation to the exchange-correlation self- 
energy at the Fermi level in Dyson's quasiparticle equa- 
tion. 

Fig. [T] illustrates how a partially filled f-band lying at 
the Ep, which is superposed on a broader band, will be 
mapped in the excitation process. If the relaxation en- 
ergy is zero or a constant, then the spectroscopic 
image of the Kohn-Sham density of states (DOS) in 
Fig. [IJa) obtained, for example, via photoemission and 
inverse photoemission processes, will be an undistorted 
copy of the DOS [l^. In general, however, E^"'^ will 
be non-zero and differ between localized and itinerant 

bands, and between occupied and empty states. There- 

( f) 

fore, it is useful to introduce the notation Ei- {P) and 
Ei^*\p) for the relaxation energies of the f and itiner- 
ant bands for the occupied states, which could be probed 
via photoemission. Similarly, Er^\l) and Er^\l) de- 
note the corresponding relaxation energies for the unfilled 
states, which could be accessed via an inverse photoemis- 
sion process. Fig. [IJb) shows that when the filled portion 
of the Kohn-Sham DOS is mapped, f-bands are shifted 



by Er {P) and itinerant bands by Ef!" (P), so that the 
f-bands no longer appear to be at the Ep- Fig. [TJc) 

shows what happens when the unoccupied portion of the 

( f) 

DOS is excited. Now the f-bands move by Ey ' (/) and 
the itinerant bands by Er^\l). The net result is that 
filled and unfilled portions of the DOS, as seen by com- 
paring panels (b) and (c), will in general appear to be 
separated in energy in their spectroscopic images. The 
preceding effects arise purely from the way the excitation 
processes play out and refiect differences in the screening 
of the added hole or electron in various orbitals of the 
unperturbed system. 

Concerning technical details, we note first that the 
ASCF scheme works well for localized core states with 
pure ^-character Q. Here, we apply the scheme to a lo- 
calized RE f-band in NCCO, where our band calculations 
based on the local spin density approximation (LSDA) 
^1, 24] and the Linear Muffin-Tin Orbital (LMTO) 
method [H-Illl, put the Nd-f band right at Ep- LSDA 
yields the correct Cu02-plane Fermi surface and reason- 
able band dispersions in doped LSCO Also, LSDA 



correctly predicts both metallic and ferromagnetic phases 
in manganites [2^ . That LSDA fails to describe the anti- 
ferromagnetic (AFM) phase within the Cu02 plane in 
undoped cuprates is not relevant here since our focus is 
on the RE physics. Notably, however, LSDA calculations 
for LSCO using lower linearization energies describe the 
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FIG. 1: (Color online) Schematic illustration of how the 
Kohn-Sham spectrum is mapped in excitation processes in 
the presence of relaxation effects, (a) Kohn-Sham DOS in 
which the filled portion is shaded green and the empty por- 
tion is shaded pink. The filled f-band is shown by a solid blue 
line and the itinerant band by a solid red line. The unfilled 
portions of these bands are marked by broken lines of the 
same color. Ef is the Kohn-Sham Fermi energy, (b) Image 
of the filled portion of the DOS where f and itinerant bands 
are shifted by Ei^\P) and Ei^^\P), respectively, in exciting 
an electron, (c) Same as panel (b), except that here the un- 
filled portion of the DOS is excited and the f and itinerant 
bands undergo relaxation shifts of Er^\l) and Ei^'^\l). 



undoped as well as the doped system reasonably [29], in- 
dicating that only small corrections are needed to the 
LSDA to bring the paramagnetic and the antiferromag- 
netic states into the correct order. For undoped NCCO, 
LSDA does predict a metallic RE AFM ground state with 
RE order in line with experimental findings [3) . A RE fer- 
romagnetic metallic solution is very close in energy, and 
it is more convenient to use for calculating excitation en- 
ergies. Accordingly, we will study excitations from the 
fictitious ferromagnetic ground state for a; = 0, which are 
very similar to the excitations of the metallic compound 
with x^Q [30l-[3^. Our basis set includes f < 3 for Nd, 
£ < 2 for Cu and O, and £ < 1 for empty spheres, two 
of which are inserted per formula unit in the most open 
part of the structure. The lattice constant, ag, is 7.45 
a.u., and the Wigner-Seitz radii are (in units of op): 0.44, 
0.33, 0.31 and 0.272 for Nd, Cu, O and empty spheres, 
respectively. The number of irreducible k-points used is 
80 for a cell containing four formula units. The com- 
puted moment on each Nd is 3.15 while the moments 
on Cu and O atoms are less than 0.1 ^b- Interestingly, 
by imposing the observed AFM order on Cu through a 
staggered field, the magnetic order on the RE does not 
change, implying weak interaction between Cu and RE 
sites. Turning to excited state simulations, we consider 
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first the process in which an occupied electron is excited 
leaving behind a hole in the system. We model this hole 
as a hole in the local density of states (LDOS) on a par- 
ticular site (e.g. Nd, Cu or O). Specifically, the hole is 
created by removing electrons from the LDOS over an 
energy window [Ec, Ep], where Ec < Ep is a, cut-off en- 
ergy defined such that a total of one electron is removed. 
The total charge for the excited state at the site t' then 
is 

Nt,iimtAE,r)dE 
-Y. I ' Nt,AE)Rl,i{E,r)dE + ^, (4) 

where Nt i{E) is the LDOS and Ri i{E,r) is the radial 
wave function at the site t with angular momentum 
The last term in Eq. |4] imposes charge neutrality within 
the simulation cell of volume f2. This form is appropriate 
for XPS, since at large excitation energies the wave func- 
tion of the electron is free-electron like with a constant 
charge density [s^. At lower photon energies (e.g. for 
UV photoemission [s^l), 1/^^ could be replaced by the 
proper density for the excited state. Calculations involv- 
ing core levels show good convergence already for small 
unit cells Q. The reason is that the charge screening 
occurs near the excited site. We have embedded one ex- 
cited site in a fairly large cell of 28 sites, and verified that 
both charges and local DOS on sites far from the excita- 
tion approach those of the ground state for all cases. The 
constraint of Eq. |4]is repeated for each self-consistent it- 
eration until the total energy is converged to an accuracy 
of about 0.01 eV. The excitation energy e„ is calculated 
by a change in total energy between the ground state and 
the excited state as in Eq. [T] 



III. RESULTS 

Before discussing our results for the relaxation ener- 
gies with reference to Table 1, we emphasize that in our 
modeling the hole does not involve a single energy level 
of the solid, but rather a group of states around the mean 
energy e(P) [sBj. The values of e(P) are given in the sec- 
ond column of the Table. e{P) is quite small (0.04 eV) for 
the majority (up spin) Nd band, dominated by f-levels, so 
that this hole state is modeled reasonably in our scheme 
[36j . Values of e(P) for other states considered (minor- 
ity spin Nd, Cu, and apical and planar O atoms) range 
from 1—4 eV, and thus the hole in these cases involves 
a substantial mixture of states around the mean value. 

Focusing on the relaxation energies E^''\P) of the oc- 
cupied states given in the 3rd column of the Table, we 
see that the up spin Nd state, which is almost exclusively 
of f-character undergoes a relaxation shift of —4.45 eV, 
so that in the excitation spectrum this state will appear 
at an excitation energy of —0.04 — 4.45 — —4.49 eV, i.e. 



at a much lower energy than in the Kohn-Sham spec- 
trum. In contrast, an excitation from the down spin Nd 
band, which contains almost no f-electron (but involves 
itinerant d-electrons) , with e = —3.63 eV will appear at 
—0.68 eV, i.e. at a much higher energy than for the up 
spin band. Along these lines, Cu excitations experience 
a relaxation correction of about —2 eV, while the correc- 
tions for planar and non-planar 0-atoms [0(1) and 0(2)] 
are both about -1-0.6 eV. As expected, the more locahzed 
levels generally suffer larger relaxation effects. In partic- 
ular, the majority spin Nd bands of f-character, which 
are most strongly localized, display the largest correc- 
tion. The Cu-d band is also relatively localized and un- 
dergoes substantial hole screening. But Cu and O states 
are hybridized so that actual excitations from the itin- 
erant CuO bands will be an average over the Cu and 
O contributions in the Table. [This is not the case for 
Nd-f bands, which hardly hybridize with other valence 
states.] Putting all this together, we estimate that an 
excitation from the majority Nd-f will be detected ^ 2.5 
eV lower than from other bands (i.e. the value of 4.45 
eV in Table 1 is reduced by about 2 eV as the Fermi 
level follows the CuO hybrid band), giving the appear- 
ance of a gap between the Nd-f states and the Fermi level, 
as shown schematically in Fig.[ljb). The aforementioned 
computed downward shift of about 2.5 eV is in reasonable 
accord with photoemission measurements, which find a 
broad peak centered at about — 3 eV below Ep for the 4f 
level |37| . 

TABLE L Relaxation energies in NCCO for excitations from 
majority- and minority-spin Nd [Nd (t) and Nd(4,)], Cu, and 
planar and apical O sites [0(1) and 0(2)]. Ei"\P) is the 
relaxation energy of photoemission with the hole at an average 
energy of e(P) and an excited electron above Ep. For the 
inverse process, the relaxation energy and the average energy 
of the added electron in the final state are Ei"\l) and e(/), 
respectively. All energies are in eV. 



Band 


e(P) Ei."\P) 


Ei"\l) 


Nd (t) 


-0.04 


-4.45 


0.04 


4.81 


Nd (;) 


-3.63 


2.95 


1.78 


4.38 


Cu 


-1.03 


-2.09 


4.54 


2.46 


0(1) 


-1.81 


0.58 


5.48 


3.40 


0(2) 


-2.96 


0.65 


5.91 


2.99 



Relaxation effects in the inverse excitation process are 
considered in the last two columns of Table 1. In this 
case, the final state involves an extra electron with an 
average energy of e(/) and the corresponding relaxation 
energy is Er"\l). Computation of the relaxation energy 
follows along the lines described above for the case of a 
final state with an extra hole, except that here we add 
an electron near the Ep in the empty LDOS on various 
atomic sites, and Ec > Ep. In order to ensure charge 
neutrality, Eq. |3]is now modified so that the charge den- 
sity term is subtracted (instead of being added) on 
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the right hand side of Eq. Turning to the resuhs for 
the majority Nd (f mainly) states in Table 1, we find 

our key result: There is a large positive shift Ei"'\l) of 
4.81 eV, implying that Nd-f(t) bands will appear to lie 
at energies well above Ep and above the itinerant CuO 
bands. In contrast to the case of photoemission, the re- 
laxation energies of Nd-t and Nd-| are very similar. This 
is because both up- and down-spin f-states exist in the 
unoccupied Nd DOS, while the occupied Nd DOS has 
almost no down-spin f-electrons. 

The Gd is a prototypical member of the rare earths 
which is a suitable benchmark for testing excitations of 
4f electron systems in X-ray photoemission (XPS) and 
bremsstrahlung isochromat (BIS) In our ground 

state calculation for the FM configuration of Gd8Cu40i6, 
the majority Gd-f band is completely filled and at about 
4.5 eV below Ep while the minority Gd-f is empty and at 
about 0.5 eV above Ep. When we apply our corrections 
the estimated XPS peak is at about —7 eV whereas the 
BIS peak is near 4 eV, which compare fairly well with the 
corresponding experimental values in pure Gd ^] , about 
—8 and 4.5 eV respectively. 

IV. DISCUSSION 

It is interesting to note that in the case of a Hubbard 
band strong on-site Coulomb replusion U between elec- 
trons of opposite spin splits the band into a lower por- 
tion lying ^ U/2 below the Ep, and an upper portion 
lying ~ [7/2 above the Ep. The gap between the up- 
per and lower Hubbard bands so created is the result of 
strong electronic correlations. In contrast, in the present 
LSDA calculation the partly filled f-band is located at 
Ep, but screening effects in the excitation process make 
the filled portion appear well below the Ep and the un- 
filled portion appear well above the Ep. The effective 
splitting between these two portions from Table 1 (first 
row), (/) - E^r"^ (P), is - 9 eV for Nd. The itinerant 
bands in Table 1 are seen to display smaller splittings. 
Apparently, the relaxation energy can be sizable even for 
itinerant valence electrons (39j . Experimental determi- 
nations of effective U for elemental Nd are in the range 
6-7 eV [H, Constrained LSDA calculations have been 
performed to estimate Hubbard U parameters fiol], and 
have been used to justify opening of Mott gaps in the 
cuprates. Similar methods have been applied for calcu- 
lation of U parameters for d-band impurities via relax- 
ation energies of embedded atoms ^41j . Our computa- 
tional scheme is however tailored for treating the process 
of electron or hole excitation for a specific pair of ini- 
tial and final states. Therefore, our relaxation energies 
cannot be described in terms of a single atomic param- 
eter, although we would expect our relaxation energies 
to be of the order of the commonly used U parameters. 
Our approach is in the spirit of the early work of Herbst 
and Wilkins [i^l although Ref. 4^ considers excitations 



from renormalized atoms within truncated Wigner-Seitz 
spheres. We also emphasize that our approach is quite 
different from LDA-fU type calculations, since we employ 
the constrained DFT formalism and not the orbital de- 
pendent techniques invoked in LDA-I-U [11]. Finally, we 
note that our LSDA calculations do not address Mott- 
physics involving strong correlations, but demonstrate 
that final state corrections can give the appearance of 
a gap in spectroscopic data even in the absence of strong 
electron correlations. 



V. CONCLUSION 

We have presented a DFT-based scheme for obtain- 
ing relaxation energies relevant for excitation of occupied 
and empty states in various spectroscopic probes. Our 
method is particularly suitable for f-bands and demon- 
strates clearly that the screening of excitations in f-bands 
differs substantially from that of more delocalized states. 
The net effect is that in the excitation spectrum the filled 
f-bands appear to move below the Ep while the unfilled 
f-bands are shifted above the Ep. The most important 
conclusion of our study is that the aforementioned shifts 
are induced via screening that occurs within the exci- 
tation process, even when the majority spin f-states in 
NCCO lie at the Ep in the LSDA Kohn-Sham spec- 
trum. We expect this screening mechanism to be ap- 
plicable to f-electron systems more generally since many 
f-electron compounds contain narrow, partially filled f- 
bands at the Kohn-Sham Ep with little hybridization 
with other bands. Although we have assumed a high 
excitation energy for the generic purposes of this study, 
it will be straightforward to extend our scheme to con- 
sider lower energy excitations. For a realistic descrip- 
tion of the spectral intensities in various spectroscopies, 
one will need to take the matrix element effects into ac- 
count. Furthermore, it will be interesting to examine 
the extent to which f-electrons in the ground state can 
contribute to Fermi surface related properties. In partic- 
ular, if NCCO has f-electrons at or near the Ep, which 
probes could detect them? Compton scattering [4J| and 
positron-annihilation [45|, |4^ , which are sensitive to the 
electron momentum density of the many-body ground 
state could be promising in this connection. 
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